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1. Introduction 

Understanding the internal structure of nucleons has long been a central mission of the theo- 
retical and experimental nuclear physics research. Not only nucleons are the building blocks of the 
visible matter, but also are they the essential probes into dynamics of the strong interaction, giving 
us glimpse of how strong interaction works to bind the quarks and gluons together to form protons, 
neutrons and other hadronic matter. Of particular interests are the nucleon form factors and struc- 
ture functions, which allow us to study the size, shape and spin structure of the nucleons. Due to 
the non-perturbative nature of the strong interaction, in which the coupling constant becomes too 
strong at low energies for perturbatiion theory to be reliable, first-principle calculations of these 
quantities are best done by lattice QCD approaches. While in the past decade or so the numerical 
techniques of such calculations have been improved greatly, some results are still plagued with sys- 
tematic errors that remain to be understood, among which chiral extrapolation, finite volume and 
excited-state contamination effects have been the most actively investigated ones recently. Most 
notably, lattice results for the nucleon axial charge, gA, are consistently lower than the experimental 
value by 10%-20%. Similar behavior has also been observed in the results for the isovector Dirac 
charge radius. While some studies indicated finite volume effects to be the culprit [1, 2], others 
found that excited-state contaminations may contribute significantly to the discrepancy [3, 4]. Of 
course, almost all the lattice simulations to date have unphysically heavy pion masses, and undeter- 
mined eiTors associated with chiral extrapolations to the physical point may constituent a big part 
of the discrepancy as well. 

It is thus to our interest to perform lattice calculations in a large volume, with a sufficiently 
large source-sink separation to suppress excited-state contaminations and at pion masses as close 
to the physical point as possible. With this in mind, we did our calculations on the 2+1 flavor 
domain wall fermion gauge configurations [5] generated by the RBC and UKQCD collaborations 
with a spatial volume of about (4.7fm) 3 , and two pion masses of about 170 and 250 MeV. As we 
will discuss later, we also set the source-sink separation to be about 1.3 fm, which is among the 
largest used in nucleon three-point functions in recent lattice calculations. In these proceedings, 
we will focus on the calculation of the nucleon electromagnetic form factors, represented by the 
Dirac and Pauli form factors F^(Q 2 ) and F^iQ 2 )- Calculations of the nucleon axial charge and 
strangeness content on the same gauge ensembles are discussed in separate presentations [6, 7] at 
this conference. 

The Dirac and Pauli form factors are defined through the matrix element 

{N( J /)\J»{x)\N(p)) = jV-P>*u(p') 

where p and p' are the momenta of the incoming and outgoing nucleons, respectively, and Q 2 = 
—q 2 = — (p' — p) 2 . It is convenient to define the isovector current (x) =J^(x)—J' j [(x) =u{x)y lx u{x) — 
d{x)y il d{x), with the corresponding isovector Dirac and Pauli form factors, Ff~"(Q 2 ) and F^' 1 ' (Q 2 ) . 

2. Details of Calculations 

The calculations were performed on the 32 3 x 64 gauge configurations with 2+1 flavors of 
domain wall fermions (DWF) and the Iwasaki gauge action with the dislocation-suppressing- 
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determinant ratio (DSDR) at /? = 1 .75. The number of sites in the fifth dimension for DWF was set 
to be L s = 32, resulting in a residual mass of am res = 0.001842(7). Ensembles with two different 
light dynamical quark masses were used, with ami = 0.001 and 0.0042, corresponding to a pion 
mass of roughly 170 MeV and 250 MeV, respectively. The strange quark mass was fixed to be 
am s = 0.045, which is close to the physical strange quark mass. The lattice cutoff was determined 
to be a- 1 = 1.37(1) GeV, or a « 0.146 fm. These gauge configurations were generated jointly by 
the RBC and UKQCD collaborations, and their basic ensemble properties can be found in Ref.[5]. 
We used the standard proton interpolating operator 

Xs(x) = e abc {[u s a {x)] T Cy 5 d s b (x)) u s c (x), (2.1) 

where the script 5 denotes the quark field smearing. In our calculation, we used Gaussian-smeared 
operator for the sources, and both Gaussian-smeared and local operators for the sinks of the two- 
point correlation functions, which we write as 



C s (t-t src ,p) = JV^Tr \&4(0\Xs(x,t)x G (0,t src )\0) 



(2.2) 



The three-point functions always have Gaussian smearing at both the source and the sink. The 
current operator, of course, remains local. We define our three-point functions as the following 

C h (0 = I^Tr[^a(0|^ G (f,^)7 M (z,0Z G (0,^)|0)], (2.3) 

where J^x) = q{x)y^q{x), ^ = (1 + Ya)/2 and ^> 53 = (1 + 74)7573/2. 

Since the nucleon correlation functions decay quickly in the time direction, we chose to use 
four different source locations on a gauge configuration to take advantage of the large time extent. 
For both ensembles, we placed the sources at t src = 0, 16,32 and 48, while the spatial coordinates 
of the sources were fixed at the origin, that is, (x,y,z) S rc = (0,0,0). The calculations were done 
on gauge configurations separated by every eight molecular dynamics units. On the ami = 0.001 
ensemble, 103 gauge configurations were used, while on the ami = 0.0042 ensemble, 165 gauge 
configurations were analyzed, making the total number of correlation functions 412 and 660 re- 
spectively. These parameters are summarized in Table 1 . 



ami 
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L 3 
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m n [MeV] 


m K L 


am res 


# of configs. 


# of meas. 


0.001 


0.045 


32 3 


x 64 


32 


170 


4.0 


0.0018 


103 


412 


0.0042 


0.045 


32 3 


x 64 


32 


250 


5.8 


0.0018 


165 


660 



Table 1: Details of the calculations. 



For the three -point functions to have sufficient overlap with the nucleon ground state, we need 
to pay particular attention to both the nucleon sources and the source-sink separations. The optimal 
choice of the source-sink separation depends on the choice of the nucleon source. As mentioned 
above, we used the Gaussian smearing for our sources. As pointed out in Ref.[8], gauge link 
smearing can help suppress the gauge noice. Hence we also applied APE smearing to the gauge 
links in the construction of the nucleon sources. Thus our proton operator consists of smeared 
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quark fields defined as 



W,„g(.x)=[ 1 --7— £ 2y{x)-U' i {x)v{x + i)-U?{x-i)y(x-i) 
\ m G 1=1,2,3 



2 



(2.4) 



where U\{x) is the APE smeared gauge link, obtained by applying N = 25 steps of the following 
transformation, 

r 4 

Ufa) = {\-c)Ufa) + - £ [U v (x)Ufa + $)Ufa+fi) (2.5) 

" v =l,v^ 

+ tft(jE-v)^(jc-v)I/ v (jc-v + j&)] . (2.6) 

We set a = 6.0, «g = 70 and c = 0.4, which we found give a nucleon operator that has a good 
overlap with the ground state [9]. 

The source-sink separation, t sn k — t src , was set to be 9 lattice units, amounting to about 1.3 fm 
in physical units. As we will see in Section 3, this choice gives us nice plateaus to determine the 
form factors for the nucleon ground state. We fixed the sink momentum to be p = (0,0,0), while 
the momenta for the operator, q = ^-n, takes all the possible values with \n\ 2 < 4. 

3. Preliminary Results 

We follow the steps as described in [10] to obtain form factors from the nucleon three -point and 
two-point functions. We first average over the four source locations before performing the standard 
jackknife analysis. No further blocking has been done yet, thus it is possible that the errors are 
underestimated. As we have only included connected contractions in our calculations, we will 
only present results for isovector quantities, in which the disconnected-diagram contributions are 
canceled. 

Figs. 1 and 2 show the quality of the plateaus for the isovector Dirac and Pauli form factors, 
respectively. For F\{Q 2 ,t), the plateaus settle in after the first two time slices for all Q 2 values, 
while for ^(Q 2 ,?), there are some fluctuations in the middle of the time slices for am; = 0.001. 
For simplicity, we choose the fit range to be t = [2,7] for all quantities. With the large statistical 
eiTors for the ami = 0.001 ensemble, the systematic errors that may result from the choice of fit 
ranges should be subleading compared to the statistical errors quoted. In the future when more 
statistics are available, the choice of the fit ranges will be re-examined. The results, labeled as 
ID32, for Fi(Q 2 ) and F 2 {Q 2 ) are shown in Fig. 3. Previous results [10] from the calculation on 
the 24 3 x 64 lattices with a lattice cutoff of cT l « 1.73 GeV are also shown for comparison (124). 
The pion mass dependence for F\ (Q 2 ) is very mild for all the ensembles, and the Q 2 dependence is 
not as steep as the Kelly parametrization of the experimental data [11], shown as the dashed lines 
in Fig. 3. The ID32 results for F2(Q 2 ), however, seem to agree with the experiment at smaller Q 2 
values, which translates into promising trend for the Pauli radii that we will discuss next. 

We also determine the root-mean-squared (r.m.s) radii from the Q 2 dependence of the form 
factors. The mean-squared radii are defined as 

2 6 dFj{Q 2 ) 

tf) = - W) ^^\&^i = h2. (3-1) 
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Figure 1: The plateaus for the isovector Dirac form factor, Ff "(Q 2 ) 
0.001 (left) and am, = 0.0042 (right). 
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Figure 2: The plateaus for the isovector Pauli form factor, FJ? "(Q 2 ), from the two ensembles with ami = 
0.001 (left) and am, = 0.0042 (right). 

We fit F\ (Q 2 ) and F2(Q 2 ) to the empirical dipole form of Fi(Q 2 ) = — fr° m which we get 

( r f) = jiji- F° r (2 2 )» as the form factors are normalized to 1 at Q 2 = 0, we only allow M\ to be 
the free parameter in the fit, while for ^(G 2 ). both ^(0) and M2 are free parameters. The results 
for the r.m.s. radii, along with results from an earlier, 124, calculation, are shown in Fig. 4. For 
(rj) 1 / 2 (Fig. 4, left), all the results seem to fall on a straight line with respect to rn 2 , and can be 
linearly extrapolated to a value which is 20%-30% lower than the experiment 1 [13]. This linear 
behavior has been observed in other lattice calculations too, as summarized in [10]. Recently, the 
authors in [4] found that the lack of curvature of the results for (r 2 ) l l 2 may be due to excited-state 
contaminations. Whether this is the case for our calculation still awaits further investigations. 

'A measurement using muonic hydrogen found a much smaller charge radius [12], which was not included in the 
PDG average. 
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Figure 3: Q 2 dependence of the isovector Dirac (left) and Pauli (right) form factors, Ff "(Q 2 ) and 
F2~"(Q 2 )- The results from the calculation discussed in these proceedings are labeled as "ID32". Results 
labeled as "124" are from Ref. [10]. 

The results for (r 2 ) 1 ^ 2 (Fig. 4, right) are more interesting. The lightest point from the ID32 
calculation and that from the 124 calculation both dip down to be away from the experiment. If 
these two points are not considered, then the results seem to approach the experiment on a smooth 
curve. One possible explanation is that the two lightest points, with m K L = 4 and 4.5, respectively, 
may suffer from large finite volume effects. Of course, there are other systematics that may cause 
this behavior, and further investigations are needed to draw a definitive conclusion. 
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Figure 4: Pion mass dependence of the isovector Dirac (left) and Pauli (right) r.m.s radii, 
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The results from the calculation discussed in these proceedings are labeled as "ID32". Results 



labeled as "124" are from Ref. [10]. 



4. Summary and Outlook 

We have reported a calculation of the nucleon isovector electromagnetic form factors on the 
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2+1 flavor domain wall fermion gauge configurations with a lattice volume of 32 3 x 64 and a lattice 
cutoff of a~ l « 1.37 GeV. Two pion masses, 170 MeV and 250 MeV, have been studied. We have 
presented some preliminary results for the isovector Dirac and Pauli form factors, as well as the 
corresponding r.m.s radii. We continue to see the linear trend in the results for the r.m.s Dirac 
radius, which, if extrapolated linearly, are about 20%-30% below the experiment. The results for 
the r.m.s Pauli radius show promising signs of approaching the experiment as the pion mass is 
decreased, with the exception of two points with the smallest m K L ( from 4 to 4.5), which may 
be an indication of finite volume effects. To understand these issues conclusively, we will need to 
increase our statistics, and perform thorough studies of excited-state contaminations, finite volume 
effects and other systematic errors. We are currently doubling the statistics for both ensembles, and 
are also applying error-reduction techniques [14] to further increase the statistical accuracy of our 
results. 

Acknowledgments 

Part of this work used the Extreme Science and Engineering Discovery Environment 
(XSEDE), which is supported by National Science Foundation grant number OCT 1053575. We 
also thank the RIKEN Integrated Cluster of Clusters (RICC) for the computer resources used in 
the calculations. M.L. is partly supported by National Science Foundation under Grant No. NSF 
PHY1 1-00905. The gauge configurations used in these calculations were generated jointly by the 
RBC and UKQCD Collaborations. The author thanks Yasumichi Aoki, Tom Blum, Chris Dawson, 
Taku Izubuchi, Chulwoo Jung, Shigemi Ohta, Shoichi Sasaki, Eigo Shintani and Takeshi 
Yamazaki for many useful discussions and important contributions. 

References 

[1] T. Yamazaki etal, Phys.Rev.Lett. 100, 171602 (2008), arXiv:0801.4016. 

[2] N. Hall, A. Thomas, R. Young, and J. Zanotti, (2012), arXiv:1205.1608. 

[3] S. Capitani et al, Phys.Rev. D86, 074502 (2012), arXiv:1205.0180. 

[4] J. Green etal, (2012), arXiv: 1209.1687. 

[5] RBC and UKQCD, R. Arthur et al, (2012), arXiv: 1208.4412. 

[6] S. Ohta, PoS LATTICE2012, 171 (2012). 

[7] C. Jung, PoS LATTICE2012, 164 (2012). 

[8] S. N. Syritsyn et al, Physical Review D 81, 034507 (2009), arXiv:0907.4194. 

[9] S. Ohta, PoS LATTICE2010, 152 (2010), arXiv:1011.1388. 
[10] T. Yamazaki etal, Phys. Rev. D79, 114505 (2009), arXiv:0904.2039. 
[11] J. J. Kelly, Phys. Rev. C70, 068202 (2004). 
[12] R. Pohl et al, Nature 466, 213 (2010). 
[13] J. Beringer et al, Phys.Rev. D86, 010001 (2012). 
[14] T. Blum, T. Izubuchi, and E. Shintani, (2012), arXiv: 1208.4349. 



7 



